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DATE / Aq^Mia S. Broitman 

^Reg. No. 38,006 

IN THE UNITED STATES PATENT AND TRADEMARK OFFICE 

In re Application of: 

Ian D. Kimber et al . : Group Art Unit: Unassigned 

Serial No.: NEW APPLICATION : Examiner: Unassigned 

International Application Number: PCT/GBOO/02434 

International Filing Date : June 22, 2000 
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Title: CORRELATION ANALYSIS 
IN THE PHASE DOMAIN 

Box PATENT APPLICATION 

Assistant Commissioner for Patents 
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PRELIMINARY AMENDMENT 

Sir : 

Prior to computation of the filing fee in the national 
stage of the above -referenced international application, please 
enter the following amendments: 
IN THE CLAIMS: 

Please amend claims 3, 4, 6, 8, 9, 10, 12, 15-16, 18,20- 
22, and 24, as follows. 
CLAIMS : 
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3. (Amended) A method according to claim 1 [or 2] , 
wherein the signals are audio signals. 

4 . (Amended) A method according to [any preceding] 
claim 1, wherein the at least one filter includes a first 
filter for suppressing frequencies which do not exhibit a 
sufficient degree of coherence. 

6 . (Amended) A method according to [any preceding] 
claim 1, wherein the at least one filter includes a second 
filter for identifying regions in the frequency spectrum of a 
cross correlation function likely to exhibit a correlated 
phase between adjacent frequencies in its Fourier Transform. 

8. (Amended) A method according to claim 6 [or 7] , 
[including] further comprising calculating the time delay 
between the common signal in the input signals by tracking the 
phase difference between the input signals as a function of 
frequency using the second filter. 

9. (Amended) A method according to claim 6 [to 8] , 
[including] further comprising calculating variations in the 
time delay between the common signal in the input signals as a 
function of frequency using the second filter. 

10. (Amended) A method according to [any preceding] 
claim 6, [including] further comprising using a third filter 
to remove frequencies which do not have sufficient amplitude. 

12 . (Amended) A method according to [any preceding] 
claim 1, wherein the at least one filter includes a fourth 
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filter for compensating the input signals for dispersion 
effects . 

15. (Amended) An apparatus [Apparatus] according to 
claim 13 [or 14], wherein the signals are audio signals. 

16. (Amended) An apparatus [Apparatus] according to [any 
of claims] claim 13 [to 15] , wherein the at least one filter 
includes a first filter for suppressing frequencies which do 
not exhibit a sufficient degree of coherence. 

18. (Amended) An apparatus according to [any of claims] 
claim 13 [to 17] , wherein the at least one filter includes a 
second filter for identifying regions in the frequency 

spectrum of a cross correlation function likely to exhibit a 
correlated phase between adjacent frequencies in its Fourier 
Transform. 

20. (Amended) An apparatus according to claim 18 [or 
19] , including calculating the time delay between the common 
signal in the input signals by traclcing the phase difference 
between the input signals as a function of frequency using the 
second filter. 

21. (Amended) An apparatus according to [any of claims] 
claim 18 [to 20] , including calculating variations in the time 
delay between the common signal in the input signals as a 
function of frequency using the second filter. 

22 . (Amended) An apparatus according to [any of claims] 
claim 13 [to 21] , including a third filter to remove 
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frequencies which do not have sufficient amplitude. 

24. (Amended) An apparatus according to [any of claims] 
claim 13 [to 23] , wherein the at least one filter includes a 
fourth filter for compensating the input signals for 
dispersion effects. 



This Preliminary Amendment is filed together with a 
Transmittal Letter to the United States Designated Office 
(DO/US) Concerning a Filing under 35 U.S. C. § 371 (entry into 
the U.S. national stage) . The above- identified amendments are 
made to eliminate multiple dependent claims and to limit U.S. 
national filing fees and to improve form. No new matter has 
been introduced. 

A clean copy of the claims is attached, as an addendum, 
hereto. 



REMARKS 



Respectfully submitted. 




Rfeg. No. 38,006 
Ostrager Chong & Flaherty LLP 
82 5 Third Avenue, 3 0*"^ Floor 
New York, NY 10022-7519 
Tel. No.: 212-826-6565 
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ADDENDUM 

SUBMISSION OF CLEAN CLAIMS 
PURSUANT TO 37 CFR § 1.121 

In compliance with 37 CFR § 1.121, the Applicants hereby 

submit a "clean" copy of the claims now pending in this 

application as follows: 

1. A method for detecting and locating a common signal 
within two input signals using correlation based techniques, 
comprising providing at least one filter by analyzing the phase 
of the input signals in the frequency domain; filtering the 
input signals in the frequency domain using said at least one 
filter; and performing cross correlation of the filtered 
signals . 

2. A method for detecting and locating leaks in a fluid 
carrying pipe using correlation based techniques, comprising: 
detecting two input signals from the fluid carrying pipe; 
analyzing the phase of the input signals in the frequency 
domain to provide at least one filter; filtering the input 
signals in the frequency domain using the at least one filter; 
and perfo3rming cross correlation of the filtered signals. 

3. A method according to claim 1, wherein the signals 
are audio signals. 

4. A method according to claim 1, wherein the at least 
one filter includes a first filter for suppressing frequencies 
which do not exhibit a sufficient degree of coherence. 

5. A method according to claim 4, wherein the first 
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filter is constructed using a method comprising: selecting at 
least one section from each of the two input signals ; 
calculating the Fourier Transform for each section; 
calculating the average vector sum of the phase difference 
between the two input signals for each of the plurality of 
frequencies; and calculating the magnitude of the vector sum 
for each frequency. 

6. A method according to claim 1, wherein the at least 
one filter includes a second filter for identifying regions in 
the frequency spectrum of a cross correlation function likely 
to exhibit a correlated phase between adjacent frequencies in 
its Fourier Transform. 

7. A method according to claim 6, wherein the second 
filter is constructed using a method comprising: selecting at 
least one section from each of the two input signals; 
calculating the Fourier Transform for each section; 
calculating the average vector sum of the phase difference 
between the two input signals for each of the plurality of 
frequencies; and calculating the magnitude of the vector sum 
for each frequency. 

8. A method according to claim 6, further comprising 
calculating the time delay between the common signal in the 
input signals by tracking the phase difference between the 
input signals as a function of frequency using the second 
filter. 

9. A method according to claim 6, further comprising 
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calculating variations in the time delay between the common 
signal in the input signals as a function of frequency using 
the second filter. 

10. A method according to claim 6, further comprising 
using a third filter to remove frequencies which do not have 
sufficient amplitude . 

11. A method according to claim 10, wherein the third 
filter is constructed using a method comprising: applying a 
digital threshold to the product of the spectra of the two 
input signals. 

12. A method according to claim 1, wherein the at least 
one filter includes a fourth filter for compensating the input 
signals for dispersion effects. 

13. Apparatus for detecting and locating a common signal 
within two input signals using correlation based techniques, 
comprising a computer including: means for providing at least 
one filter by analyzing the phase of the input signals I the 
frequency domain; means for filtering the input signals in the 
frequency domain using said at least one filter; and means for 
performing cross correlation of the filtered signals. 

14. Apparatus for detecting and locating leaks in a 
fluid carrying pipe using correlation based techniques, 
comprising: detectors for detecting two input signals from the 
fluid carrying pipe; a computer including means for analyzing 
the phase of the input signals in the frequency domain to 
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provide at least one filter; means for filtering the input 
signals in the frequency domain using the at least one filter; 
and means for performing cross correlation of the filtered 
signals . 

15. Apparatus according to claim 13, wherein the signals 
are audio signals. 

16. Apparatus according to claim 13, wherein the at 
least one filter includes a first filter for suppressing 
frequencies which do not exhibit a sufficient degree of 
coherence . 

17. An apparatus according to claim 16, wherein the 
first filter is constructed using a method comprising: 
selecting at least one section from each of the two input 
signals; calculating the Fourier Transform for each section; 
calculating the average vector sum of the phase difference 
between the two input signals for each of a plurality of 
frequencies; and calculating the magnitude of the vector sum 
for each frequency . 

18. An apparatus according to claim 13, wherein the at 
least one filter includes a second filter for identifying 
regions in the frequency spectrum of a cross correlation 
function likely to exhibit a correlated phase between adjacent 
frequencies in its Fourier Transform. 

19- An apparatus according to claim 18, wherein the 
second filter is constructed using a method comprising: 
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selecting at. least one section from each of the two input 
signals; calculating the Fourier Transform for each section; 
calculating the average vector sum of the phase difference 
between the two input signals for each of a plurality of 
frequencies; and calculating the magnitude of the vector sum 
for each frequency. 

20. An apparatus according to claim 18, including 
calculating the time delay between the common signal in the 
input signals by tracking the phase difference between the 
input signals as a function of frequency using the second 
filter. 

21. An apparatus according to claim 18, including 
calculating variations in the time delay between the common 
signal in the input signals as a function of frequency using 
the second filter. 

22 . An apparatus according to claim 13 , including a 
third filter to remove frequencies which do not have 
sufficient amplitude . 

23. An apparatus according to claim 22, wherein the 
third filter is constructed using a method comprising: 
applying a digital threshold to the product of the spectra of 
the two input signals. 

24. An apparatus according to claim 13, wherein the at 
least one filter includes a fourth filter for compensating the 
input signals for dispersion effects. 
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CORRELATION ANALYSIS IN THE PHASE DOMAIN. 

This invention relates to the detection of a common signal within two (or more) • | 

input sign-ils using correlation based techniques, and is more specifically, | 

although not excliisively, concerned with thft detectioTi and location of leaks in ' ^ 

water mains using correlation based techniques. \ | 

There arc many applications in which signals are analysed using the technique 

of crosscorrclation of two data streams. This technique is specifically useful / • ' v 

when analysing or comparing two or more composite time sequential signals • , * ' ''\ 

within whi(;h it is believed tliat tlicrc is a coinnaon signal. ' ; . |i 

One area in which crosscorrclation is used is in the analysis of waicr leak |; 
soimds conducted along water mains, specifically locadng Icalui un waici- ' [i 

■ . I':. 

mains using audio sensors. The technique involves the detection and locad^ J_ 

a common leak sound in the signals from two listening devices. It is anticipated . . / / . } 

. f 

that the sound of a leak heard remotely at one location will contain significant ;j 

5 

similarities to the souiid of the same leak heard remotely at a second location. ? 
The two audio signals are theretbre likely to exhibit a peak in their 
crosscorrclation function if there is a common source present, the position of 
the peak giving the time delay when the signals are similar. The peaks in the 
crosscoixelation function can therefore be used to deduce the location of 
common audio signals via the velocity of the soxmd in the pipe and the physical 
surroundings of the noise source. If the nature and persistency of the soxind is 
indicative uf leaks Xlien there is a Ixigh probability that a leak has been found 
and located- 
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Tlac analysis techniques described herein both foe tlic piiur art methods and the 
present invention aie also applicable in iiimiy oLhcr technical areas. These 
include the following; the analysis of odd mechanical noise in equipment; the 
assessment of quality of objects by the noise that they make while travelling 
along a production brie; the assessment of the. performance of wheeled vehicles 
on road or railway tracks; the analysis of road or rail surt'aces by the analysis of 
the sound of the passage of vehicles either on a mobile platform or at a fixed 
location; and the measurement of closely similar recordings (the analysis of 
their true relationship from a high quality digital copy to a different 
performance of the same song that sounds very similar). 

Coii elation processing tccluiiques are fiequeiilly pciToiincd a^inv^ die following 
process; calculate die Fourier ti'ansforni of die data seiii rcprciicndng the two 
wavefurmii to be cuirclatcd. thereby transforming the data into the frequency 
domain; performing the required mtilxipUcations; and calculating the inverse 
Fourier transform to display the correlation function in the rime domain. Often 
the Fast Foinier Transfomi is used since this allows a considerable reduction in 
the processing required to generate the correlation fiinction. These techniques 
are well known and generally used. 

The general correlation processing techniques work well, ond several devices 
are available to perform this taak. However, when the signal to noise ratio is 
poor ond/or the bandwidth of the signals is restricted due to the propagation of 
the signal then the crosacotrdation function inevitably becomes degraded, with 
the peaks becoming less distinct or even hidden by the noise. Thus fillers are 
used lu obifciin liic optimum signal to noise ratio and coirclaiion accuracy. 
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Noise which may be present in a signal can be considered as two distinct 
forms: random noise and clutter. Random noise is unpredictable on a ehort 
timcscalc, but has a zero mean. It may be removed by averaging over a 
sufficiently long period. Clutter is a function of the detected signal such as, for 
example, echoes and resonance, which cannot be removed by averaging. 

The signal to noise ratio may be improved by fihcring the input signals to 
reduce the bandwidth. However, as the bandwidth is reduced the width of the 
correlation peak tends to increase. This makes ii difficult to determine the time 
delay when the signals are similar, and thereby makes it more difficult to 
deteraiine the location of the common signal. It is therefore imrortant tn filter 
the signal to exclude only those bands which do not contain the useful signal. 
However, the prior art methods of selecting the bands are imprecise. 

For example, in the case of water leak detection the audio noise signal picked 
up &om the pipe will be a combination of the leak noise and other external or 
waterbomo soimds- In order to unprove the signal to noise ratio of the input 
signals the opemtor of a water leak detector may analyse the spectrum of the 
input signals to determine where signal power is significant. Filters arc then set 
to remove die frequencies in which the signal power is not significant. This 
mcLhod uatially provide some improvenients, but does not guarantee a band 
with good correlation. The opcrdior will gcncraJly need to uac a combination of 
experience and mal and error to improve the shape of the correlation fiincrion. 
This process can be time consuming, often taking tens of minutes at each leak 
sire. 
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Tlie present invention seeks to mitigate these disadvantages, and to provide an 
improved mcdiod of detection and location of a signal using correlation based 
techniques. 

According lo a first aspect of ihe invention, there is provided a mcthud for 
detecting and locating a common signal within two input signals using 
correlation based techniques, comprising providing at least one filter by 
analysing the phase of the input signals in the frequency domain; filtering the 
input signals in the frequency domain using said at least one filter, and 
perfomung crosscorrelation of the input signals* 

According to a second aspect of the invention, there is provided a method for 
detecting and locating leaks in a fluid carrying pipe such as, for example, a 
water main using correlation based techniques, comprising detecting tvvo input 
signals firom the fluid carrying pipe- analysing the phase of the input signals in 
tlic frequency domain to provide at lea^t one filter; filtering the input signals in 
the frequency domain using die at least one filter, and performing 
crosscorrelation of die inpui signals. 

According lo a third aspect of the invention, there is provided apparatus 
for detecting and locating a common signal withia two input signals using 
correlation based techniques, comprising a computer including: means for 
providing at least one filter by analysing the phase of the input signals in the 
frequency domain; means for filtering the input signals in the frequency 
domain using said at least one filter; and means for perfomung crosscorrelation 
of the filtered signals. 
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ActurdiuR to a fourth aspect of the invention, there is provided apparatus for 
deteciijiy, diid locating leaks in a fluid carrying pipe using correlation based 
techniques, compri5iinj=;: detectors for detecting two input signala from the fluid 
canning pipe; a compuLcr including means for analysing the phase of the input 
STgTi.?ils in The frequency domain to provide at least unc filter; mcaiis for 
filtering the input signals in the frequency domain using the at least ouc filter; 
and means for performing crosscorrelarion of the filtered signals^ 

The signals maybe audio signals. 

The at least one filter may include a first filter for suppressing frequencies 
which do not exhibit a sufficient degree of coherence. The first filter maybe 
cousiructed using a method comprising: selecting at least one section from each 
of the two input signals; calctilating the Fourier Transform for each section; 
calculating tlie average vector sum of the phase difference between the two 
input signals for each of a plurality of frequencies; and calculating the 
magnitude of the vector sum for each frequency. 

The at least one filter may include a second filter for identifyiag regions in a 
frequency specimm of a cross correlation function likely to exhibit a correlated 
phase between adjacent frequencies in its fotirier transfomi- The second filter 
may be constructed using a method comprising: selecting at least one section 
from each of the two input signals; calculating the tbuiier transform for each 
section; calculating the average vector sum of the phase difterence between the 
two input signals for each of a plurality of frequencies: and calculating the 
magnitude of the vector sum for each frequency. 
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The ihne delay between the common signal in the input signals may be 
calculated by tracking the phase difference bctvir-ccn the input signals ns a 
functional frequency using Llic second filter. The variation in the time delay 
between the conmion signal in the inpuL isiuuals as a functional frequency may 
be calfAilatftd using a second filter. 

There may be a third filter for rftmovjng frequencies which do not have 
sufficient amplitude. The third filter may be constructed using a method 
comprising applying a digital threshold to the product of the spectra of the two 
input signals. 

The at least one filter may include a fourth filter for compensating the input 
signals foi dispersion effects. 

A specific embodiment according to the invention is provided for the detection 
of leaks in water mains. Tlic apparatus may comprise a novel correlator 
provided on a sUmdard field engineer's Personal Computer, The novel 
correlator may include a suite uf signal processing algorithms and a suitable 
Graphical User Interface. Like cuirent correlators the novel correlator would 
receive its inputs from two n:ansducers. Input signals are processed by the 
computer using signal processing techniques including those described herein 
to detect and locate leaks. 

As explained hereinbefore, the characteristics of a correlated signal are that it 
has similar features in its waveform (although one may be delayed with respect 
to each other and the whole buned in noise and interference). The 
CTOSSCorrelation of r%vo input signals will only produce a significant peak if the 
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frequency bands which niakc up ihost; input signals have a stable phasic 
relationship with respect to each other, i.e. the input signals have a narrow hand 
cohercTice. Filtering of the rwo input signals to remove or suppress the 
frequencies which do not exhibit a sufficient degree of coherence will therefore 
enhance the peak in the crosscorrelation fimction. Frequencies which do not 
exhibit sufficient coherence need not be included in the final result because 
they will only be contributing noise. The filter protile for the signal may 
therefore be determined automatically. 

The novel coixclator includes a phase confidence filter, which is conatructed to 
suppress ihc fiequeacics iu the input signals which do not exhibit a sufficient 
degree of coheienctj. The pha^c i;unfidcucc fiher effectively indicates the 
estimation of how well a parricular frequency will uorxclalc, and is the result of 
analysing the phase stability between the input signal^;. One pus:>iblc m.cihod 
for constructing the phase confidence filter would be JorsiL lo calculate a pha^c 
confidence function using the following steps: to calculate the complex Fourier 
transform for each of a number of sections of the two input signals; to calculate 
the average vector sum of the phase difference between the two input signals 
for each frequency; and then to calculate the magnitude of the vector sum and 
normalise the result. Thus a function is provided which has a value of 1 .0 for a 
given frequency which shows perfect coherence, and tends to 1/x/n if the 
phase differences for a given frequency are random, where n is the ntimber of 
samples taken. The phase confidence function may then b© used as an optimal 
ficqucncy weighting function to construct a filter which will include only those 
frcqucueies wliich contribute usefully to the crosscorrelation of the input 
signals. 
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The phase confidence filter is used in combination with an automatic band 
filter. The automatic band filter is constructed based purely on the ampUmde of 
tlic input signals. The filter is designed to select frequencies of large amplitude, 
wliich arc likely to be prominent in the crosscorrelation, and by reducing the 
bandwidth to hence impiovc die signal to noise ratio. This filter may be 
constructed by simply applying a digital tlircshold to the product of the spectra 
nf the input signals. The threshold may be configuiablc by tlic user. The phase 
confidence filler and the automatic band fiher arc apphcd iicparately to tlic 
data. 

The use of the automatic band tilter addresses one practical disadvantage of the 
phase confidence filter. This is that some of the noise m the mput signals can 
have highly correlated noise as a result of imperfect digitisation, this effect, 
seen tisually at high frequencies, can produce a substantial phase confidence 
measiorement. The automatic band filter suppresses frequencies where the 
signal to noise ratio is poor and this artifact is most likely to be seen. 

Thc effect of using the phase confidence filter in combination with the 
automatic band filter is shown in Figs, I to 4. Figs. 1 and 2 show a first input 
signal 1 and a second input signal 2 respectively, both in the frequency domain. 
The crosscorrclalion function 7 for tlie unfiltercd first and second signals is 
shown in Fig. 3. The phase confidence filter 3 and the automatic band filter 4 
are constnicied for the first and second signals and applied to the signals in the 
frequency domain to produce a filtered first signal 3 and a filtered second 
signal h. The crosscorrelation function 8 for the filtered first and second signals 
is shown m Fig. ^. The imprnvement benveen the crosscorrelation functions 7 
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and S can clearly be sccu. Tlierc is improved clarity, with reduced background 
jftVRl and less ambiguous pcsik locauoii. 

If the conunnn features in die input signals cuutaiu some degree of time 
dependency, then the phase relationship will not be stable. The use of the phase 
confidence filter calculated using the method above may blui tlic results. An 
alternative method for calculating the pbase confidence filter would be to 
consider the phases for each individual section of the input signals. The process 
could then adapt if the phaee appeared to be stable to letigthen the averafieby 
induding further sections to improve the estinjate. 

This roctliod could also be adapted to track the signals if they were time 
dependtrnL. Consider, for example, a signal which is very aoisy but has a 
particular value winch may be either constant or changing overtime. If this 
value is consliuit then it is possible to average the signal to produce the value 
accurately. Huwevei, if the value is changing then a compromise must be made 
between the number of samples over which the signal is averaged, and the rate 
at which the value is v aryiug. 

If to a first approximation this variation is linear it is relatively easy to produce 
a root me.an square best srraight line fit thai continually updates as die data 
comes in. If the variation is non-linear then some means must be found of 
releasing the linear constraint but .still producing a good curve fit. If the 
calculation process also produces two alternative straight line fits, one based on 
the first half of the data, and the other based on the second half of the data (they 
could be updated every other data point to ensure an even number of pninrs) if 
the two lines show coherent deviation above some threshold a track deviation 
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can be irtitiated. Thereby, an adapiive signal iiiicking process may be 
developed. 

The novel correlator also inclndes a filler based on a consideraxion of the phase 
relationshipE between adjacent treqiiencifts. Relationships exist between 
adjacent frequencies in the different phase analyses, ihat is, as well as having 
a reasonably stable value, the phase differences may exhibit a steady 
progression of phase over the range of frequencies which are coherent The 
source of tills result can be visuahscd by considering a single correlation peak, 
a peak in ihc liiiie domain is created by frequencies adding up in phase in the 
frequency domain. Cunsidei firstly a centre correlation peak (often the result of 
an instrumental anefact or inLcrfcriii.^ signal) wliich is located in the centre of 
the correlation range and which corresponds to correlation where there is zero 
rel?itive delay between two input signals. In this case tlie frequencies have the 
same or very similar phase, and therefore the phasic docs not change with 
frequency. A centre correlation peak may be removed by coastructing a filter to 
remove frequencies having the same phase. Note, however, that exact centre 
correiaticn in the leak, although Unlikely, is not impossible. 

If a peak in the cross correlation function represents a delay between the input 
signals then frequencies in the cross correlation fimcrion of those signals are 
shifted in phase by an amount diat is proportional to the hme delay and the. 
frequency. Frequencies evenly spaced in a Fourier Transfomi wiU exhibit a 
linear chemgc of phase between adjacent correlated phases. 



Cuus>ideiing the phases of sequential frequencies can have the advantage when 
ihc prupa^iHion baindwidth is very narrow resulting in a broad correlation peak 



18-DEC-2001 16=11 PROM EDUfARD EUfiMS EPR .«J Olt^lRi^litS Q^^^^l: 

WO 00/79425 PCT/GBOO/02434 

11 

because the phace progression can be ''unwrapped" and a straight hne fitted to 
it. 

The cAisience of a second major peak in the correlation result will add a second 
linear progression lu the piiase pattciii. The two phase progressions will add as 
vectors with the reiiult thai iii> ihc rclalivc phases of the two linear progressions 
move in and out of phase the main peak will be advanced and retarded in phase 
givmg a ripple in the phase results. 

Alternatively, the case where the velocity of the sound frequency dependent 
may be utihsed. In this case the phase progression for adjacent frequencies In a 
single correlation peak will not be linear. This non-lmearity may be tracked. 
This will allow either the compensation of results for dispersion where the 
coixclation is spread by its effects or the identification and measurement of the 
dispersion effects using known sound sources. 
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CT.ATMS. 

1 . A methoci for delecting and locating a tuuniion si.^al within two input 
signals using coiXftlatioTi based techniques, compiiiiing pruvidiiig at least 
one filter by analysing thft phase of the input signals in the frequently 
domain; filtering the input signals in Xhet frequency domain using said at 
least one filter; and performing crosscorrelation of the filtered signals. 

2- A method for detecting and locating leaks in a fluid carrying pipe using 
coirelation based techniques, comprising: detecting two input signals trom 
the fluid carrying pipe; analysing the phase of the input signals in the 
frequency douiain to provide at least one filter; filtering the input signals in 
the frequency domain using the at least one filter; and performing 
crosscoixclaLiuu of tlic filtered signals. 

3. A method according to claim 1 ox 2, wherein the signals arc audio signals. 

4. A method according to any preceding claim, wherein the at least one filter 
include?; a first filter for suppressing frequencies which do not cjdiibit a 
sutficient degree of coherence. 

5. A method according to claim 4, wherein the first filter is consmicted usir^g a 
method comprising: selecting at least one section from each of the two input 
signals; calculating the Fourier Transform for each section; calciilating the 
average vector sum of the phase difference between the two mput signals 
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for each of a plurality of frequencies: and calculating the magnitude of the • ■ 

vector sum for each frequency. 

6. A meihod dcuuxdiug to any preceding claim, wherein the at least one filter ' • 
includes a second filter fui identifying regions in the frequency spectrum of • 

a cro<;scorrelation function likely to exliibit a corrclaied phase between 
adjacent frequencies in its Fourier Trauiifurra. 

7. A method according to clam wherein the second filter is coni>trucicd 
using a method comprising: selecting at least one section from each of the 
two input signals; calculating the Founer Iransfotm for each section; 
calculating the average vector sum of the phase difference between the two 
input signals for each of a plurality of frequencies; and calculating the 

magnitude of the vector sum for each frequency. • ^ > 

8. A method according to claim 6 or 7, including calculating the time delay . J' 
between the common signal in the inout signals by tracking the phase 

difference between the input signals as a function of frequency using the 
second filter. 

9. A meihod according lo any of chxiimi 6 to 8, including calculating variations 
in the time delay between the common signal m the input signals as a 
fimction of frequency using the second filter. 



10, A method according to any preceding claim, including using a third filter to 
remove frequencies which do nor have .Sufficient amplitude. 
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1 1 . A method according to c[3irr) 1 0, wherein the third filler is consirucLcd 
using a method comprising: applying a digital threshold to the product of 
the spectra of the two input signals. 

12. A method according to any preceding claim, wherei_n the at least one filter 
includes a fourth filter for compensating the input signals for dispersion 
effccLb. 

13. Apparams for dcieciiug and locating a common signal within two input 
signals using correlation b-dscd lecliniqucs, comprising a computer 
including: means for providing at Icust one filter by analysing the phase of 
The input signals in the frequency domain; meanii for filtering the input 
signals in the frequency domain using said at least one filter; and means for 
peiforming crosscorrelaiion of the filtered signals. 

14. Apparatus for detecting and locating leaks in a fluid cauying pipe using 
correlation basied techniques, comprising: detectors for detecting two input 
signals from the fluid carrying pipe; a computer including means for 
analysing the phase of the input signals m the frequency domain to provide > 
at least one tilter; means for filtenng the input signals in the frequency 
domain using the at least one filter; and means for performing 
crosscorrelation of the filtered signals. 

15. Apparatus according to claim 13 or 14, wherein the signals are audio 
signals. 
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1 6. Apparani.. according to any of cl^iim. 13 to 1 5, wherein xhc at least one 
filter inciiides a first filter for suppressing^ fiequcncics which do no: exhibit 
a sufficient degree of coherence. 

1 7 . Apparatus according to claim 1 wherein the first filter is cunbUacted using 
a method comprising: selecting at least one section from each of the two 
input signals; calculating the Fourier Transform for each section; 
calculating the average vector sum of the phase difference between the two 
Input 5ignab for each of a plurality of frequencies; and calculating the^ 
magnitude of the vector sum for each frequency. 

18. An apparatus ducurdia^ lo any of claims 13 to 17, wherein the at least one 
filter includes a second filter fur identifying regions in the frequency 
spectnim of a crosscorrclaLiun fuiit;tiuii likely to exhibit a correlated phase 
between adjacent ficqucneicb iu its Fourier Transform, 

19. An apparams according to claim IS. wherein ihc second filter is coiisu-uctcd 
U5ing a method comprising: selecting at least one section from each of ihc 
two input f^ignals; calculating the Fourier Transform for each section; 
calcuiatmg the average vector sum of the phase difference between the two 
input signals for each of a plurality of frequencies; and calculating the 
magnitude of the vector sum for each trequency. 

20. ^Vd apparatus according to claim 18 or 19, including calculating the time 
delay between the common signal in the inout signals by tracking the phase 
difference between the input signals as a function of frequency using the 
second filter. 
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21. An appRTatii5; according to any uf ddiiiis 18 to 20, including calculating 
variation? in the time delay between the commuii signal m the mput signals 
as a function of frequency using the second filter. 

22. Apparatus according to any of claims 13 to 2 1 , including a third filter to 
remove frequencies which do not have sufficient amplitude, 

23- Apparatus according to claim 22, wherein the third filter is constructed 
using a method comprising: applying a digital threshold to the product of 
Qic spectia of tlic t>?vo input signals. 



24. AppcLralus aucoidiiig to any of claims 13 to 23, wherein the at least one 
filler includes a £ouirh filter for compensating the input signals for 
dispcTbion effects. 
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ABSTRACT 

Method and apparatus for detecting and locating leaks in a fluid carrying pipe using 
correlation based techniques is described. Two input signals are detected from the fluid carrying 
pipe, the phase of the in put signals in the frequency domain is analyzed o provide at least one 
filter, the input signals are filtered and cross-correlation is performed on the filtered signals. 
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